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Abstract 

The Jacobs-Rebbi equation arises in many contexts where vortical 
motion in two-dimensional ideal media is investigated. Alternatively, 
it can be derived in the Abelian Higgs field theory. It is considered 
non-integrable and numerical solutions have been found, consisting of 
localised, robust vortices. We show in this work that the equation is 
integrable and provide the Lax pair. The exact solution is obtained 
in terms of Riemann theta functions. 

1 Introduction 

Studying the interaction energy of vortices of the Ginsburg-Landau model of 
superconductivity, Jacobs and Rebbi [Tj have derived a nonlinear equation 
from which the scalar function (order parameter) can be calculated. The 
derivation is based on the similarity of this model with the Abelian Higgs 
theory |2j, [H], j3j In a special case, corresponding to a particular choice of 
parameters, the nonlinear equation takes a simple form, a nonlinear elliptic 
differential equation in two spatial dimensions. The same problem is treated 
by Dunne [S] in a general context of gauge theories with Maxwell and/or 
Chern-Simons terms in the Lagrangean density. The same Abelian Higgs field 
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theory is discussed and the special case mentioned above is identified as the 
minimum of the action functional (saturation of the Bogomolnyi inequality), 
realized by fields obeying a simpler set of equations. This special case is 
called self-duality, with reference to the equality of the differential two-form 
(the curvature of the fibre bundle) with its Hodge dual in the geometrical 
setting of the field-theoretical content. 

In fluids and plasmas coherent motion and in particular vortices are ubiq- 
uituous They can appear even in turbulent states. In two-dimensions, 
apart from the dynamical equations derived from the conservation laws an 
alternative model has been proposed, consisting of the motion of discrete, 
point-like vortices in plane interacting via a potential 7\. It has been shown 
that this model can be mapped onto a field-theoretical model whose struc- 
ture is very similar with the non-Abelian gauge- Higgs field theory e 8j. The 
self-dual state of this field at stationarity is precisely the asymptotic state 
of the ideal fluid, which effectively provides an analytic derivation of the 
sm/i-Poisson equation describing the fluid streamfunction. For the station- 
ary states attained at very large time by the ideal ion instability in plasma 
(described by Hasegawa-Mima), the model of discrete vortices interacting in 
plane introduces a short-range potential (then a massive photon of the gauge 
field in the field theoretical model). 

There are two differences between the field theoretical models developed 
starting from the plasma problems and the model from which Jacobs-Rebbi 
equation is derived. First, the absolute value of the scalar Higgs field is 
constant at large distances (on a circle of very large radius); this means 
that the vorticity should be constant at large distance, in the plasma case. 
Second, the model must be Abelian, which is less than we would need for 
treating the case of the Euler fluid (the sm/i-Poisson equation). However, 
there are physical situations where the boundary conditions for the vorticity 
are compatible with the formulation of the Jacobs-Rebbi model. And, it 
is known that Abelian models can provide description of plasma problems 
(guiding centre particles, for example) leading to the Liouville equation. The 
effect of these differences still needs to be investigated, but in any case, the 
exact determination of the solution can only be a useful instrument. 

It is usual to consider that the Jacobs-Rebbi equation cannot be solved 
analytically and in consequence numerical solutions have been provided. 

We show in this paper that the Jacobs-Rebbi equation is exactly inte- 
grable. We consider the integrability on periodic domains and provide the 
Lax pair. We follow the standard algebraic-geometric method of integration 
and generate explicit solutions in terms of Riemann theta functions. 
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2 Derivation of the Jacobs-Rebbi equation 



2.1 Derivation in the context of Ginsburg-Landau the- 
ory 

The Ginsburg-Landau theory is the framework in which the Jacobs-Rebbi 
equation has been derived, since the original aim was the investigation of 
the energy of interaction between two vortices in superconducting Helium. 
The interest for vortical structures of Ginsburg-Landau comes also from the 
observation that the theory of a gauge field coupled to a scalar field (Abelian 
Higgs field) can in some cases exhibit also coherent vortical structures. We 
include in this Section the derivation according by Jacobs and Rebbi pp. In 
fluid physics other approaches can be developed to arrive at similar forms of 
the equation. 

The free energy of the Ginsburg-Landau theory and the potential energy 
of the Abelian-Higgs theory is 

E = j d 3 x j~%F ij 

+^|(d,-*e^|Vc 4 (|0| 2 - C 2) 2 } 

where is a complex scalar field, A{ is the Abelian gauge potential and 

Fij = d{Aj djAj 
The minimum of the energy is attained for 

101 = c ^ 

The variables are rescaled as 

1 _ 

X% X{ 

c e 
Ai = c Ai 

(j) = C (j) 

and the energy becoms 



where 




The model is restricted to the case where all field functions does not 
depend on the third coordinate and 

A 3 = 

In plane the coordinates are expressed by complex variables 

Z = X\ + %X 2 
Z = X\ — IX2 

and the differential operators can be defined 
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Analoguous combinations are used for the two remaining potential compo- 
nents 

A = \{M-iA 2 ) 
A = + %A 2 ^j 

The energy per unit length along the coordinate axis x 3 is 

e 



£ = — [ dzdz\2\dA-dA\ 



+ 



(d - 1A) 



+ 



(d - iA) 



2 A 2 



In the following the tilda will be omitted. 

The equations of motion that are derived from the Lagrangean 
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(d - iA) (d -iA)(j)+(d- iA) (d-iA)(P- — 
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4ddA - Ad 2 A 

—i(j)d(f) + i 
-2A(fp 
= 

For a general value of A one can replace particular forms for the two 
functions <p and A and obtain differential equations. 
For the value 

A = 1 

the situation is different since one can obtain a lower bound for the energy. 
The following integration by parts is done 



/■ 
/■ 



dzdz [{d - iA) <p(d + iA) 0] 
dzdz [(d-iI)^(d + iA)(t> 

-i (dA - dA) #] 

Then one obtains the new expression for the energy 



£ = I J dzdz^\(d-iA) 0| 2 + 
J dzdz (dA - dA) 
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{dA-dA)--(\<p\ 2 -l) 



Taking into account the boundary conditions and asking for the absolute 
minimum to be attained the terms in the curly braket must be taken zero 
and we obtain the equations 

(d - iA) = 

dA-dA+^(\(f)\ 2 -1) =0 
This equations can be further transformed 

A = idip 
A = —idip 

then the first equation reduces to 

(9 — dip) (j) = exp (ip) d [exp (— ip) <fi] = 
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which means that we can introduce an analytic function 



/ = exp (--0) (p 

Inserting 

(z, z) = exp [-tp (z, z)} f (z) 
in the second equation we obtain 

9^ = ^ [exp (2^) //-l] 

By a new substitution 

i> = X - 2 ln iff) 

the equation becomes 

9dx = g [exp (2x) - 1] 
with the condition the x goes to zero at infinity. 



3 General procedure for obtaining solutions 
of the Jacobs-Rebbi equation 

The procedure is similar to those developed for the sine-Gordon equation [5| 
and for sinh-Poisson equation Just as in the general case of nonlinear 
differential equations which are exctly integrable by the algebraic-geometric 
procedure, we start from a configuration which is specified initially. By con- 
trast with other equations (for example KdV, etc) where the variables are 
space and time and the unknown function is given at t — 0, here the coor- 
dinates are both spatial. Then the conditions to be specified are boundary 
conditions, for example taken on the lines x = and y = 0. 

Consider that a certain flow configuration is specified (for example from 
experimental measurements) and the boundary conditions are specified in 
the form of two functions 

u 0x (y) &ndu 0y (x) 

for x G [0,L] and y G [0, L\. We assume that L is much smaller than 
the radius of the circle on which the asymptotic value of the vorticity is 
given. The "initial" values of the unknown function is introduced in the 
Lax operator eigenvalue problem. Solving this problem we identify a set 
of eigenvalues (the Lax operator spectrum, see ^T]) and the corresponding 
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eigenf unctions with periodicity properties (Bloch functions). It is a general 
situation that in the spectrum there is a subset of eigenvalues for which the 
two eigenf unctions are identical. These eigenvalues are called non- degenerate 
and the subset is called main spectrum. 

Using the main spectrum one can construct the hyperelliptic Riemann 
surface associated with the Wronskian of the eigenfunctions. For a two by two 
Lax operator, (i.e. hyperelliptic Riemann surface) the monodromy problem 
is simple. 

One has to define on this surface the dual homological sets: cycles and 
differential one-forms. Then the period matrices can be calculated. 

Using the inverse of the v4-period matrix one can generate the variables 
(the phases) appearing in the arguments of the Riemann theta function. 

Finally, one can calculate the solution at any point (x, y) and can repre- 
sent it graphically on a space domain. 



3.1 The spectral problem for the Jacobi-Rebbi equa- 
tion on periodic domain 

From a detailed consideration of Lax pairs found by Forest and McLaughlin 
for the sine-Gordon equation we obtain the following Lax equations. The 
first is 



A- _ ff. —nA 1 ( du i -du 

V r l 8x 4 ydy dx 

This set of equations is considered on a periodic domain along the x axis. It 
is of second differential order and has two independent solutions periodic on 
x, which we note and 0_. We chose them to correspond to the following 
initial conditions at X — Xq 




<p + (x ,x ,p) = 
(p-(x ,x ,p) 



A nv other solution of the system, corresponding to the following condition 
taken at x = x 

(j)(x = Xo,p)= ( q 

is a linear combination of these two basis functions 

(X, Xq, p) = ( ^ j^' ^| J = P0+ (X, X ,p) + Q0_ (X, Xq, p) 
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We consider the second set of equations 

Ui % p ~ VP ~ dy ~~ 4 ( 

A 2 



du I ■ du 
dy ^ l dx 




This is a system of two differential equations on the periodic domain along 
the y axis, having two independent solutions. These two functions must 
actually be identical to the previously defined functions, <f> + and <p_ since 
finally we can only accept solutions of both sets of equations, on x and on y. 
On the y direction we take the initial conditions 



<l>+(yo,yo,p) = 
(j)-(yo,y ,p) = 



Any other solution of the system, corresponding to the following condition 
taken at y = y 

( P' 

^(y = 2/o, P) = ( q, 
is a linear combination of these two basis functions 

The fact that we write only the x or the y notation is only derived from the 
context of the first or second system. Actually the pair of functions (<p, ip) 
depend on (x, y) on a two-dimensional periodic domain and are independent 
solutions of the two systems of equations. We note 

.du\ 



16^ 
„• 9 , 



- Vp 



dx 
A 2 

— + iw 

dy ^ 4 W 



y/P 



_ (du 

\dy dx J 

-i-k + > 

exp (u) - y/p 




= 



or 



dx 



-W(j>i + 



Vp 



-dx- + 



16^ 



exp (u) - y/p 
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exp (u) - y/p 02 = 
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16VP 
^ + 4^ + 



16^ 

v 7 ^ J V>i 



.9^2 i 



A 2 



exp (u) - v /p 



1^2 
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dijj 2 i , 



16VP 



exp («) - v /p 
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According to standard procedures we define squared eigenfunctions, as 
combinations of the components of the two independent solutions and -0. 



/ = 

9 = 
h = 



2 

Ml 
-^202 



(-0102 + 01^2) 



and calculate the derivatives at x and at y. 
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dg 
dy 

dh 
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16VP 
: exp (u) + ^/p 



Vp 



160? 



exp (u) 



w 

f ~2 9 



w 



+ VP~)f + ^h 



Vp- 



exp (u) 



u ■( r- A2 

h + l I vT~ 



16,/pJ 



iw 

— o + 2? 
2 y 



A 5 



iw , 

— ft + 2? 



I^Vp 
.( A 2 



exp (u) - ^/p 



/ 



V16VP 



- Vp) f 



(i) 



(2) 



(3) 



Using the squared eigenfunction it is possible to construct the constant of 
motion 

C = f-gh (4) 

with the properties 

dC_ 

dx 
dC 

dy 

Then C depends only on the eigenvalue p 

C = C(p) 

It can be shown that the Wronskian of two solutions of the systems of equa- 
tions 



W = det 



4>i 4>i 

4>2 1P2 

= - 02^1 



can be expressed in terms of this constant of motion by the relation 

W 2 

The fact that we have a formal expression for the Wronskian in terms of 
a function of only the eigenvalue p, allows us to discuss the problem of the 
existence of two independent solutions to the systems of equations. There 
will be independent solutions everywhere on the complex p plane except at 
the points where the Wronskian vanishes. The set of points on the complex p 
plane where the Wronskian vanishes (and there is only one solution) is called 
main spectrum of the scattering problem. 

We will write the squared Wronskian (i.e. C) as a polynomial of the 
variable p thus formally introducing the points of the main spectrum, Pi,i = 
1,27V. 

2N 

"4^ = 11^ -») ^ 

i=i 

Since there is a relation between the squared eigenfunctions and the Wron- 
skian, we will introduce analoguous expressions as polynomial in the variable 
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p 

N 



f = ^X>< (.) 



N 



9 = ^9kP l 

k=0 

N 



fc=0 



where the coefficients are functions of (x,y). 

We dispose of differential equations relating these functions, Eqs.(|2|) and 
(J3J), and we will insert the polynomial expansion and find relations between 
the coefficients. 



k=0 

,2 N 



+ \y^/§ Yl + exp (u) h k ] p h 



It results 



1 A 2 

= -— [go + h exp (u)] (7) 

dfi 1 . N 1 A 2 r , . , 

-^T = 2 (fo + ^o) + 2 16 ^ + 1 6XP ^ 

dfk 1 1 A 2 

~q~ = 2 ^ 9h ~ 1 + hk ~^ + 2 16 ^ 9k + hk 6XP ^ ' k = 2 ' '"' N ^ 

o = ^ (g N + h N ) (9) 
Using now the definition of the Wronskian and the polynomial expressions 

C = f-gh (10) 



2N 

} I (P - PA; 



fc=l 
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In order the invariant quantity C (p) (or the Wronskian) to be a polynomial 
in p, there should be no source of singularity in Eq. fjlOj) and this means that 
the function / (p) must have 

/o = 

as we have already taken in (J2J). 

The coefficient of the zero-degree of p in C (p) is 



C(j>=0) = -g h 



Ea. ffTUj) gives 



27V 



C(0) = (-l) 2N l[p k = -g h 



(11) 



(12) 



k=l 



We now introduce the zeros jk (x,y) ,k = 1, N of the function g (x, y;p) 
(we suppress the arguments (x,y)) 



N 



9 (p) = n & - ^) 



(13) 



fc=l 



from which it results 



N 



9(0) = (-If II 



lk = go 



(14) 



k=l 



From the equations JEJ), (??) and ((T3J) we obtain 
exp (u) : 



^0^0 



1 ^ A^ rr^ 
-1) Hfe=i7fe 



l-r2JV 



or 



u = In 



ll*=i7fc 

T-T27V 
llfc=lPfc 



(15) 



(16) 



Eq. (fTB^) shows that if we have the main spectrum and if we could calculate 
the zeros of the squared eigenfunction g, we could find the solution to the 
nonlinear equation. The ensemble of zeros of the squared eigenfunction g is 
called auxiliary spectrum. 
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3.2 The equations for the auxiliary spectrum 

To find the differential equations obeyed by 7^ we start from the Eqs.(|2} for 



9g 2 

dx 



exp (u) + y/p 



w 

f ~2 9 



and its y version 



dg iw 

0y- = ~Y 9 + 2i 



exp (u) - y/p 



f 



and calculate all terms at p = 7^, a zero of g, using Eq.(JT3J) 



dg_ 

dx 



P=7k 



A 2 



exp (tt) + v /t^ 



f(lk) 



(17) 



(19) 



The value of / (7^) is determined from the expression of C (p) , Eqs.Q and 
(fTUj). after inserting p = j k 



2N 



n(^-^)=[/(7*)f 



In Eq. ()19|) we will also replace exp (u) from Eq. ()15j) . Then 



dx 



IL=i H 



nffi (7* - Pi) 



1/2 



nfti (7* ~ t») 



In a similar way we have 



dy 



N 



n ~ 



2i 



i=i 
A 2 



16V7fc 



exp (u) - 



or 



^7fc 

ay 



-2i 



y-rN 
lL=l7l 



16 nf=iw 



1/2 



n£i (7fc - 70 



(20) 



(21) 
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3.3 Checking the formulas as solutions 



As Ting, Chen, Lee ^H] have shown for the case of the sm/i-Poisson equation, 
it may be useful to try to find out if the formulas determined above, Eqs.(j20J) 
and (f2*T]) may already be taken as solution, for a set of which is not yet 
determined. The procedure consists of replacing the expression (fTBj) in the 
initial equation, perform the derivatives of the functions 7^ (x, y) appearing 
in this expression and taking into account the equations of motion, Eqs. ()20|) 
and (HQ- 

The following change of variables makes the calculation easier 



y -»• V 

and the initial equation becomes 

d 2 u 



x — > x = ix + y 
—ix + y 



X 2 

— [exp (it) - 1] 



dx'dy' 

The equations of motions are translated in the new variables 



dx' 



2i 



dy' 
X 2 



y-rN 
IL=1 H 



16 VTk uZpi 



+ V 7fe 



nSi (7* - pi) 



1/2 



Uti (lh - 7/) 

l=t=k 



d-f k d~f k 



dx 1 
-2i 



dy' 
X 2 



il/=i 7z 



n£i (7*: - pi) 



1/2 



16 VTk uZpi 

Adding and substracting these equations we obtain 



Y[i=x {ik - n) 

l=f=k 



dx' 



2iV7fc 



Ul=i (ik - pi) 



1/2 



n£i dk - 11) 

l=f=k 



dlk 
dy' 



X 2 (nili7z) Ul=Alk-Pi) 



1/2 



8 vfk uZpi n=i(7,-70 

l=f=k 



(22) 
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This last expression can be written, using Eqs. (jl5|) and 



dy> 



A 2 i ( , d lk 



The conversion formula is 



exp {u) 



T~f N 



Ui=iPi 

and can be used to obtain the derivatives of u 



— = 2exp(- W 



N 



Ui=iPi , =1 

Il^Li (7i - Pm) 



N 



z=i 



y- 1 #7; 



Ilm=l (7/ - 7m) 



(23) 



(24) 




Using Eqs.(J22J)> (jZHI and the conversion equation (J2"4*)) this expression is 
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rewritten 



d 2 u 



dy'dx' 16 / rf 2 \dx' J ^(li-Pm) 

1 / A 2 \ , , 1 fdji\^y 1 /1<9tz 1 #7 



~SB)(- 



A_ 
16 



exp (it, 



1 








\ " 

m=l 


1 


H - 




dli\ 




dx' J 



2i V 16/ GXP U/ Jji \dx' J ^ii-im \li dx' <y m dx' 



We finally obtain four terms in the expression 



d 2 u A 2 

- — exp (u) (Ti - 2T 2 + 2T 3 - T 3 ) (25) 



dy'dx' 16 

JV 1 / Q . \ 2 2Af 



£?7? V^V ^(ll-Pm) 



N 1 / Q N 2 AT , 

07M 1 



n/ 2 ( A™/ ) 



(=1 ' N ' m=l 



m^l 

N N 



\- \- 1 1 d'Jl d'Jm 

{^~[ll- 1m lllm 9x' dx' 



m^l 

S 



y- , 

' 7, 3 \ 



«97« 



This expression must be compared with 

d 2 u A 2 



dy'dx' 8 



[exp (it) - 1] (26) 



It can be verified that, for an arbitrary set of pt,k = 1,...2N, and a 
set of functions ji (x, y) , I — 1, N verifying the differential equations (j2~U|l 
and (j2H) (or, equivalently, Eqs.fl22) and (EH ) the two 

expressions f)25j) and 

()26|) are identical. This means that the initial nonlinear equation is verified 
if u(x,y) is given by the expression (|16jl. The verification can be done by 
summing the residuues in a formal expression defined by integration in the 
complex plane of a function having an adequate singularity structure. We 



16 



note however that the expression can be verified also on purely algebraic 
grounds, chosing arbitrary sets {pk} and {'Jk}- From these numbers one 
calculates the derivatives appearing in the four terms of the above formula, 
without any need to solve the differential equations. The expression is ver- 
ified simply as an algebraic expression, by a symbolic software, or, for any 
particular choice the verification can be done numerically. 

We conclude that we dispose at this moment of a method to find a solution 
of the Jacobs-Rebbi equation on a periodic spatial domain. This consists of 
chosing a set of 2N arbitrary complex numbers, {pk} and solving the first 
order differential equations for {'jk} with a set of initial conditions. 

3.4 Solving the equations for the auxiliary spectrum 

To solve the differential equations for 7^ (x, y;p) , k = 1, N starting from a 
set of initial conditions is a difficult task as is apparent from the form of the 
Eqs. (J2Uj) and (|7T|). However there is a standard procedure that provides the 
analytic solution of these equations. It is based on the fundamental property 
of 7^ (x, y; p) of being defined when p maps the complex plane (of the spectral 
variable of the Lax operator) to the complex function given by the square 
root of the Wronskian. Since the later is a polynomial in p, the square root 
defines a hyperelliptic Riemann surface, i.e. a compactified double covering 
of the complex plane with cuts connecting pairs of zeros of the Wronskian. 
These are the points {pk, k = 1, 2N} of the main spectrum, plus the point 
zero and the point at infinity. The point zero appears since in formulas (j'2()j) 
and (|2T|) a factor of can be adjoined to the product of the 1 = 1,..., 2N 
differences (7^ — pi), simply by taking formally p^ = 0. Then the object 
which can be defined on the basis of the square root of the Wronskian but 
reflecting the need for the particular form in the equations of 7^ 's is 



with po = 0. The geometry of this hyperelliptic surface is important in 
finding the solution. 

Pairs of zeros pk are joined by cuts and in addition the origin is connected 
to infinity. This gives a number of iV + 1 cuts and generates a compact 
Riemann surface of genus g = N. 

On this surface there are defined two objects characterising the differential 
geometry of the curve: 

• a basis of the one dimensional cohomology group of the surface; this 
means two sets each of iV closed paths on the curve (cycles), having 



2iV 
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particular intersection properties. The two sets are noted aj, and re- 
spectively bj, j — 1, N. The intersections are 

aj o dk = 
a . ob k = 5 jk 
bj o b k = 

A typical example, for an elliptic curve g — 1 with the topology of the 
torus, consists of the two possible closed turns around the torus, the 
short way (a) and the long way (£>). 

• a basis in the ring of the one-dimensional differential forms 

^ dp 

d^h = p , v , k = 1, ...,iV 

With these two sets one calculate several quantities which are invariants 
of the Riemann surface. Essentially there are calculated integrals of the 
elements of the basis of differential forms along the cycles aj and bj. These 
are called periods and are organised in two matrices 



= / dfii = 

J a-; J a-: 



Aj = I dm = / ,i=l,N,j = l,N 



It is useful to work with the inverse of the matrix A 

C = A- 1 



Using C, the matrix of A periods is reduced at the identity matrix, and the 
matrix B becomes 

r = CB (27) 

the r- matrix, with positive imaginary part. 

Using this geometrical framework the solution of the 7^ equations can 
be obtained by oparating first a transformation from the set {7^} to a set of 
functions {<ftk} representing phases of motion along the cycles of the Riemann 
surface. This transformation effectively linearises the motion, which can be 
trivially integrated in these new variables. 

We have to define the functions of the target set, the phases {(/>&}. They 
are integrals of linear combinations of the differential one-forms along paths 
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on the Riemann surface, each starting from an initial point 7 and ending in 
the point which correspond to a function 7;. The integrand is a combination 
of the differential one-forms with coefficients from the matrix C = A~ l 



/11 j< 
, -1 '0 m=l 



(2f 



The mapping that realises the correspondence from a collection of points 
{7/, 1 = 1, N} of the hyperelliptic Riemann surface to a manifold defined by 
the collection of points {(f>k,k = 1,N} is called Abel map. The manifold 
generated by the points k = 1, N} has genus g = N (as the initial curve) 
and has the topology of a torus. It is called Jacobi torus. 

Since the upper limit in the integrals are precisely our points 7^, we can 
obtain the differential equations for If. by direct derivation of this formula 
and using the differential equations for 7*,. 



N a N 
°li \- 



dx 



A^-faA^ C krad^m (ll) 
1=1 m=l 



N 



N 



^2 C km22 



and 



m=l 



N 



* R (71) dx 



dy 



m=l 



U R (li) dy 



Replacing the derivatives from Eqs.(j20J) and (|2TJ) we have 



dx 



and 



N 



2 ^ Ck m ^ 



ll 



N-m 



1 E[n=l {ll - In) 

nj^l 



m=l 



n 



N 

n=l In 



+ Vii 



(29) 



N 



N 



dy 



2? ^ Ckm ^ 



ll 



N-m 



m=l 



tt Un=l (ll 

n=t=l 



In 



X 2 



rin=l In 



1Q Vli UT =1 Pn 



(30) 

We have to calculate separately the two terms in each of the above for- 
mulas. 

N-m-l A 2(rin=l7n, 



ll 



1=1 lln=l (ll -In) 16 rin=l Vn 

n=£l 
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N „ N-m 



i=i n«=i (7/ - in) 

Since the product of all the 7„ 's is independent of the summation index I, we 
will factorse it, as well as the product of the eigenvalues p n and the constant 

2 



A 2 (]ln=l7nJ ^ 7 iV-m-l 



pl ~ m " — 2^ 



16 Il^lPn l=1 Un=l ill ~ In) 

Tracy (for the case of Nonlinear Schrodinger Equation) and Ting, Chen, 
Lee (for sm/i-Poisson equation) (TH] adopt different procedures to calculate 
the sum. For example, one can write the Lagrange interpolation fromula for 
an arbitrary function on a set of N points {xk} 



f( X )=2^f ( X 3) 
3=1 



Y\ k =l (Xj - X k ) 

k+j 



Then one takes 
then 



f{x) =7* 



n rifc=i ( x - x k) 



, 3 Y\k=l (xj - x k ) 

If the product at the numerator is expanded one gets a polynomial of degree 
N while in the left hand side we have a polynomial of degree q. Comparing 
the coefficients of the same powers of the variable x in both sides it is obtained 

N 1 

% - (jv - 1)] = J>; 



, , nf=i ( x j - x k) 



From this we find that 



N JV-m-l 

^Twh 7 = S[N-m-l-(N-l)} 



~[ I!"-! ill - In) 



S (m) 



and 



r = S[N-m-(N-l)} 



1=1 Un=l ill ~ In) 



8 (1 — m) 
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Pi 



and the equations becomes 



A 2 



n 



N 

71=1 



In 



16 UZlPn 

p 2 = 5(l-m) 



-5 (m) 



dx 



N 



m=l 



km 



A 2 



n 



A' 



71=1 ^ n 



-5 (m) + 6(1 — m) 



dfa 
dy 



N 



2i ^ ^ Cfcm 



m=l 



a 2 ill;) i 
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-2iC kl 



t-{2N 

Un=lPr. 



-6(m)-6(l 



m) 



The equations can be trivially integrated and we obtain the (x, y) depen- 
dence of the phases 



fa 0, y) = 2C fe i (x -iy) + 4> 



kO 



(31) 



where fao are constants of integration, initial phases. 

We note from Eq.(|3ip that the motion on the Jacobi torus is entirely 
determined by the main spectrum through the topological properties of the 
hyeprelliptic Riemann surface (canonical cycles, differential forms, period 
matrices). 



3.5 The Jacobi inversion 

After the determination of the phases fa, which are points on the Jacobi 
torus, we want to be able to retrive the functions jk (x, y) of the auxiliary 
spectrum, since they are necessary for the explicit determination of the solu- 
tion u (x,y), via Eq.(|16p. This constitutes the Jacobi inversion problem and 
has been solved in connection with elliptic functions. The main instrument 
is the Riemann theta function. 

The definition of the Riemann theta function involves a vector of dimen- 
sion N (we denote it by <fi) and a N x N matrix r whose elements have the 
imaginary part positive. 

oo oo / N N N \ 

((f), r) = ^2 • ' ' ^2 eXP I ^ ^2 mk< j )k + ni m i T V m i ) 

mi=— oo m n =— oo \ k=l i=l j=l / 
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In general the G function is associated to a hyperelliptic Riemann surface of 
genus N generated for example from a two sheeted covering of the complex 
plane with 2N + 1 or 2N + 2 branch points, between which N + 1 cuts ahve 
been done. The matrix r corresponds to the matrix determined from the 
periods of the canonical differential one-forms on the canonical cycles, see 
Eq. (j27|) . The argument of the G function is the vector . The following 
periodicity properties of the G functions are useful in the inversion problem: 

1 . The translation with unity of only one component of the vector <fi leaves 
G invariant. This is actually related with the fact that the components 
of the arguments <pk are coordinates along the cycles of the g-torus, and 
so they are periodical. Using the symbol for a column vector of N 
components having only one 1 in position k and in rest, we have 



2. Adding to the argument <fi a vector consisting of one of the columns 
(say, j) of the matrix r generates a factor to the function G 



The function G with argument a vector of dimension N has N zero's. 
These roots of the Q function solves the Jacobi inversion problem. 

Certain necessary quantities must be defined. We consider again a linear 
combination of the canonical differential one-forms dfik with coefficients taken 
from the columns of the matrix C = A' 1 . These linear combinations are 
integrated on the Riemann surface along paths staring from an arbitrary 
point 7o and ending in some point of the surface, q 



These are functions of the current point on the Riemann surface, q. We 
consider the sum of the integrals of such functions along the a-cycles plus 
terms from the diagonal of r 



G(0 + e fc ,r) = G(0,r) 



Q + Tp r) = exp (27rz0j — irirjj) Q ((f), r) 




(32) 





(33) 
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Finally, it is considered the function 



C(g) = e(z/(g)-0 + D,r) 



(34) 



It is proved that the zero's of the function ( (q) are 7&, the auxiliary 
spectrum. 



3.6 Solution of the Jacobs-Rebbi equation in terms of 
Riemann ^-functions 

Any initial condition for the nonlinear Jacobs-Rebbi equation leads to a main 
spectrum, i.e. a set of complex numbers {pk,k = 1,2N}. From these we 
construct the hyperelliptic Riemann surface of genus iV and calculate the 
period matrices and the phases of the linear motion along the canonical 
cij cycles on the surface. This is purely topological and geometrical data, 
generated from the main spectrum or, equivalently, by the initial condition 
for the unknown solution u. 

On the other hand, solving the Jacobi inversion problem provides us with 
the auxiliary spectrum {7^, k — 1, N} where 7^ are functions of the phases, 
and, as such, of the variables (x, y). 

The eigenvalues of the main spectrum and the functions 7^ (x, y) of the 
auxiliary spectrum give the explicit form of the solution u(x,y) via the con- 
version formula 



u 



In 



t-iN 

Hjfc=i7fc 



21 



n 



2N 



I Pm 



(35) 



N 



2N 



2^1n 7fc - J^ln 



Pn 



k=l 



m=l 



Returning to the result of the Jacobi inversion procedure, we will try to 
express the first sum in terms of the function's zero's, i.e. in terms of the 
zero's of the function ( (q). 

In Ting, Chen and Lee [TU] it is adopted the method consisting of gen- 
erating directly the sum of the logarithms from an integral of a complex 
function. 

We have to remind that the hyperelliptic Riemann surface is a mapping 
from the complex plane of the spectral parameter p via the square root of the 
polynomial expression generated by the Wronskian. The variable q appearing 
as the upper limit of integration in Eq. ()32)1 is a point on the hyperelliptic 
Riemann surface and is the image of a point in the complex p-plane; the 
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path of integration in Eq. (j32j) is the image of a path on the same p-plane. 
We can try to introduce an intermediate object, a complex function whose 
singularities will lead us, after integration, to the sum of logarithms. This is 

with the contour of integration T being a path on the Riemann surface. This 
path must be chosen such that it circles all the zero's of £ (q), and then the 
value of the integral will be 

TV 

where qoj,j = 1,N are 

C(%) = o 

The contour F is specified after the Riemann surface is mapped back onto 
the p-plane as the normal polygon obtained from cutting along the canonical 
dj and bj cycles. Since the genus of the Riemann surface is N, the number of 
cycles is 2N and each cycle generates two edges of the polygon, with opposite 
senses. The polygon has 4N edges and it is chosen as the contour T. All the 
points qoj are somewhere inside the polygon, so what we need is a choice for 
the function / (q). It is natural to take 

f(q) = \nq 

since we want the sum of the logarithms, but this induces an additional 
singularity at q = and a cut connecting q = to oo on the p-plane. This 
cut on the p-plane is translated into two paths on the hyperelliptic Riemann 
surface and since the variable of integration on the path is £ we have to 
connect the two points into which p = is mapped with the single point on 
the surface that corresponds to p = oo. This actually separates the polygon 
T into two closed parts. The integration in Eq. (f3T)j) must be done separately 
on the two contours. A part of the integration will be done along the cut and 
in one integration the path correspond to one determination of the logarithm 
(one branch) while in the other integration the path is on the next branch of 
the logarithm. 

f°° dC f°° dC 
J+ = - / hig-^+ / (lng + 27rz)-± 
Jo s Jo s 

2m I ^ 
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[°° (In q + 2m) ^ + f° (In g + 4ttz) ^ 
Jo s Jo <■> 



/•°° rfC 

The rest of the integration is the sum over the poles of the integrand, i.e. 
the zero's of the function ( (q) 

These are two images of the compactified two sheeted covering of the 
complex plane in two Riemann hyperelliptic curves with genus 2 and respec- 
tively genus 3. They correspond with the case where the number of branch 
points in the main spectrum pk is 2g + 2 = 6 and respectively 8 , i.e. if the 
eigenvalues comes in three or four pairs. 



Figure 1: Two sheet Riemann surface of y — y/x 2 — 1. 

An alternative calculation of the same integral Eq. (j36|) is done following 
directly the path along the canonical cycles a and b. In this evaluation the 
periodicity properties of the G function are essential. When the point of 
integration q is on a b k cycle, the point that corresponds to it but attached 
to the opposite side of the cut along the cycle can be reached by a complete 
turn along the nearest a cycle. But such a change does not introduce any 
modification in the integrand, since it is exactly the operation involved in the 
first periodicity property of 0. This means that the integration along the 
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Figure 2: Branched covering of the the complex plane by the two sheet 
Riemann surface of the function y = y (z 2 + a) (z 2 + b) (z 2 + c). This would 
correspond to a main spectrum consisting of only n = 6 points. 



Figure 3: The hyperelliptic Riemann surface of genus g = 2. This is 
topologically equivalent to the two-sheet Riemann surface shown in El since 
n = 2g + 2 = 6 gives g = 2. 
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Figure 4: A hyperelliptic Riemann surface of g = 3. 



edge of T coming from a cycle b k can be paired with the integration along 
the edge coming from the opposite side of the cut along b k , without any 
change in the integrand. Since these two integrals are equal but of opposite 
sign we conclude that the edges originated from 6-cycles do not contribute to 
the integral. Analoguous consideration for the a^-cycles involve the second 
property of the function. If q is on an edge representing one side of the cut 
along the a k cycle, the point q on other side can be reached by a complete 
turn along a b cycle. This introduces the change of the integrand 

C(?) C (?) = ex P {-i^Uk ~ 27ri [u k (q) - (j) k + D k }} C, (q) 
This means 



In C (q) = -mr kk - 2%i [u k (q) - 4> k + D k ] 
+ \n((q) 



and 



d\n( 



c 

-2mdv k (q) + In C (q) 



Taking into account the conclusion reached before that the b cycles do not 
contribute to the integration and leaving aside the part coming from the 
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branch cut the integral becomes 



1 [, d( 

2m J T C, 



N r 

E 

fc=i 

E 

k=l J a k 



d((q) 



a k S 



dC(?) 



lno 

a k S 



In g 



dC(?) 



In q 



2mdv k (q) + 



d((q) 



or 



1 

2th 



(3? 



This integral is a constant since <iz/fc is a differential one-form generated from a 
linear combination of the canonical one-forms (depending only on the surface) 
and the integration is performed over closed loops a^. It does not leave any 
choice since it does not depend on any parameter. 

We have completed the calculation of the integral (pTfj) in the two ways: 
one with the polygonal dissection of the Riemann surface plus the branch 
cut (which gives the right hand side of (jrS7j) ) and one with the path on 
the surface, using the reunion of canonical cycles, obtaining the constant of 
Eq. (j38|) . It only remains to make explicit the last term in Eq. (j37|) coming 
from the branch cut integration. 



dC(q) 



In 



CM 
C(o) 



(39) 



with the relation 



CM 



e\6 + v(oo)- Dl 



e 



+ z/(0)+ / dv(q)-T> 
Jo 



In the argument, v (0) — D is a constant that can be included in the initial 
phases 4>ko (Eq. (j31j) ). The integrals of the differential forms are done along a 
path that can be completed with a circle at infinity. It results a loop can then 
be mapped onto the set of loops that surround the cuts, i. e. effectively it is 
shrinked the set of a-cycles. The integrals are then reduced at the diagonal 
entries of the A matrix which are all unity. 

1 f 1 
dv (q) = - / dv (q) = -I 

J c 
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Then Eq. fl39|) can be written 



<*C(g) 

c 



In 



0(0) 



We return now to the Eq. (jH7j) and (|3*5j) 



1 

2tu 



lng 



E ln ^- 2 / 

fc=i 
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where from we obtain 

N 



j=x - ^ k=l Ja * 

The explicit form of the solution is given by the conversion formula (|35|) 



N 



2N 



a 



2^1n 7fc -^ln 



Pr, 



k=l 



m=l 



4 In 



e (0+ii) 



6(0) 

fc=l ^ afc m=l 



Since the last line is composed of constants, 

N 2N 



k=l ^ ak m=l 



we can write the solution as 

u (x, y) = 4 ln 

with 



0(0) 

0fc (x, y) = 2C fc i (x - iy) + feo 
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Figure 5: The solution of the Jacobs-Rebbi equation. 



Fteu .4log[e(z l + 1/2,z 2 +1/2,Z3+1/2;,);e(z l ,z 2 ,z 3 ;i)] 




Figure 6: The solution clearly shows the vortical structures as expected. 
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4 Conclusion 



In conclusion we have proved that the Jacobs-Rebbi equation is exactly inte- 
grable and have provided the exact solution. We have followed the standard 
approaches developed in detail for similar equations: sine-Gordon and sinh- 
Poisson equations. 

Knowledge of the exact solution will make more accessible the investiga- 
tion of the physical applications of this equation. 
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